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Within the framework of the constituent quark model, it is shown that the single hadron frag- 
mentation function of a parton can be expressed as a convolution of shower diquark or triquark 
distribution function and quark recombination probability, if the interference between amplitudes 
of quark recombination with different momenta is neglected. The recombination probability is de- 
termined by the hadron's wavefunction in the constituent quark model. The shower diquark or 
triquark distribution functions of a fragmenting jet are defined in terms of overlapping matrices of 
constituent quarks and parton field operators. They are similar in form to dihadron or trihadron 
fragmentation functions in terms of parton operator and hadron states. Extending the formalism to 
the field theory at finite temperature, we automatically derive contributions to the effective single 
hadron fragmentation function from the recombination of shower and thermal constituent quarks. 
Such contributions involve single or diquark distribution functions which in turn can be related to 
diquark or triquark distribution functions via sum rules. We also derive QCD evolution equations for 
quark distribution functions that in turn determine the evolution of the effective jet fragmentation 
functions in a thermal medium. 
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I. INTRODUCTION 

In the study of the properties of dense matter and 
search for quark gluon plasma in high-energy heavy-ion 
collisions, jet quenching 1] or the suppression of leading 
high pt hadrons Q has become a powerful diagnostic 
tool. The theoretical concept relies on the study of par- 
ton propagation in medium and induced radiative energy 
loss @H II El- Using perturbative QCD (pQCD) and 
taking into account the intricate Landau-Pomeranchuk- 
Migdal (LPM) interference effect, one can calculate 
the radiative energy loss. It has a unique dependence 
on the path length traversed by the parton and the local 
gluon density of the medium. Such energy loss is man- 
ifested in the modification of the effective parton frag- 
mentation functions |8| which in turn leads to the sup- 
pression of the single inclusive hadron spectra. In non- 
central heavy-ion collisions, the path length dependence 
of the energy loss gives rise to the azimuthal angle depen- 
dence of the high pr single hadron spectra ||. Therefore, 
experimental measurements of the suppression of single 
hadron spectra or the modification of the fragmentation 
functions and their centrality dependence can provide im- 
portant information about the initial gluon density and 
geometry of the produced dense matter. Moreover, one 
can go beyond single hadron spectra and study the mod- 
ification of multiple hadron correlations inside jets due 
to multiple scattering of the parton in the medium |10|] . 

Data from experiments at the Relativistic Heavy-ion 
Collider (RHIC) have indeed confirmed the predicted fea- 
tures of jet quenching One has seen not only a sig- 
nificant suppression of single inclusive high pt hadron 
spectra 1 1 2L Il3| and strong azimuthal angle dependence 
[l4j , but also the suppression of high px hadrons on the 
opposite side of a triggered high px hadron in the cen- 



tral Au -f- Au collisions 0] . These are all consistent with 
the qualitative features of jet quenching due to parton 
energy loss |l6[ . The extracted initial gluon density in 
the most central Au + Au collisions at y/s = 200 GeV 
is found to be about 30 times higher than that in a cold 
Au nuclei [Til UtI IT^ | . Combined with the enormous col- 
lective expansion as measured by the elliptic flow of the 
final bulk hadrons [l9l | , current experimental data points 
to the formation of a strongly interactive quark-gluon 
plasma in central Au + Au collisions at RHIC. 

During the propagation and interaction inside a decon- 
fined hot partonic medium, a fast parton can have not 
only induced gluon radiation but also induced absorption 
of the surrounding thermal gluons. The detailed balance 
leads to an energy dependence of the net energy loss that 
is stronger than without for an intermediate energy par- 
ton poj . In principle, one can consider gluon absorption 
as a parton recombination process and it continues until 
the hadronization of the bulk partonic matter. Even- 
tually, during the hadronization, partons from the jet 
can combine with those from the medium to form the fi- 
nal hadrons. Indeed, there already exists some evidence 
for such parton recombination in the measured hadron 
spectra in heavy-ion collisions at RHIC. At intermediate 
Pt = 2 — 4 GeV/c, the suppression of baryons is signifi- 
cantly smaller than mesons, leading to a baryon to meson 
ratio larger than 1 21]. This is about a factor of 5 in- 
crease over the value in p + p collisions. On the other 
hand, the azimuthal anisotropy of the baryon spectra is 
larger than that of meson spectra. Such a flavor depen- 
dence of the nuclear modification of the hadron spectra 
and their azimuthal anisotropy is not consistent with a 
picture of pure parton energy loss followed by vacuum 
fragmentation. The most striking empirical observation 
of the flavor dependence, that could reveal the under- 
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lying hadronization mechanism, is the scaling behavior 
between the azimuthal anisotropy of baryon and meson 
spectra j^, vf{pr/2)/2 = uf (p r /3)/3. Such an obser- 
vation is inspired by a schematic model of hadron pro- 
duction by constituent quark recombination 23). Here, 
D2 («a ) ^ s the second coefficient of the Fourier trans- 
formation of the azimuthal angle distribution of mesons 
(baryons). 

Many quark recombination models [24l I25I l2l| are suc- 
cessful in describing the observed flavor dependence of 
the nuclear modification of hadron spectra at intermedi- 
ate pt- These models in general involve thermal quarks 
in the medium and employ the constituent quark model 
for hadron wavefunctions which determine the recombi- 
nation probabilities. They, however, differ in the han- 
dling of specific recombination processes. Some consid- 
ered only recombination of thermal quarks, with inherent 
correlations caused by jet quenching. Others also include 
recombination of thermal and shower quarks from the 
fragmenting jet, which dominate the hadron spectra at 
intermediate pt- They also differ in the determination of 
the constituent quark distributions from high px jets and 
there exist ambiguities in the connection between partons 
from pQCD hard processes and constituent quarks that 
form the final hadrons. 

One of the models that we will follow closely in this 
paper is by Hwa and Yang |24|. In this model, quark 
recombination processes are traced back to parton frag- 
mentation processes in vacuum. They assume that the 
initial produced hard partons will evolve into a shower 
of constituent quarks which then recombine to form the 
final hadrons in the hadronization process. The for- 
mulation of the recombination of the shower quarks of 
the partonic jets and the medium quarks in heavy-ion 
collisions is straightforward, given both the shower and 
medium quark distributions. Since the Hwa- Yang model 
is a phenomenological one, the nuclear modification of 
the shower quark distributions and their QCD evolution 
cannot be calculated. The model has to rely on fitting 
to the experimentally measured hadron spectra to ob- 
tain the corresponding nuclear modified shower quark 
distributions for each centrality of heavy-ion collisions 
and correlations between shower quarks are completely 
neglected. 

In this paper, we make a first attempt to derive the 
recombination model of jet fragmentation functions from 
a field theoretical formulation and the constituent quark 
model of hadron structure. Within the constituent quark 
model, we consider the parton fragmentation as a two 
stage process. The initial parton first evolves into a 
shower of constituent quarks that subsequently will com- 
bine with each other to form the final hadrons. Since 
constituent quarks are non-perturbative objects in QCD 
just like hadrons, the conversion of hard partons into 
showers of constituent quarks is not calculable in pQCD. 
However, we can define constituent quark distributions 
in a jet as overlapping matrices of the parton field op- 
erator and the constituent quark states, similarly as one 



has defined hadron fragmentation functions. Ignoring in- 
terferences in the process of quark recombination, we 
demonstrate that the single inclusive meson (baryon) 
fragmentation functions can be cast as a convolution of 
the diquark (triquark) distribution functions and the re- 
combination probabilities, which are determined by the 
hadrons' wavefunctions. This is similar in spirit to the 
Hwa- Yang recombination model. Given a form of the 
hadrons' wavefunctions in the constituent quark model, 
one can in principle extract constituent quark distribu- 
tion functions from measured jet fragmentation func- 
tions. We are also able to derive the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) IH H<| evolution 
equations for the diquark and triquark distribution func- 
tions, which in turn give rise to the DGLAP evolution 
equations for the hadron fragmentation functions within 
the quark recombination formalism. Such a reformula- 
tion of the jet fragmentation functions does nothing to 
simplify the description of jet hadronization. However, 
extending the formalism to the case of jet fragmentation 
in a thermal medium and working within field theory at 
finite temperature, we can automatically derive contri- 
butions from recombination between shower and thermal 
quarks in addition to soft hadron production from recom- 
bination of thermal quarks and leading hadrons from re- 
combination of shower quarks. The shower and thermal 
quark recombination involves single (diquark) quark dis- 
tribution functions which can be obtained from diquark 
(triquark) distributions through sum rules. Therefore, 
one can consistently describe three different processes 
within this formalism. One can also consistently take 
into account parton energy loss and detailed balance ef- 
fect for jet fragmentation inside a thermal medium. 

For convenience in this paper, we will simply refer to 
the constituent quarks as quarks and the initial current 
quarks and gluons as partons. The quarks from jet frag- 
mentation will be denoted as shower quarks in contrast 
to quarks in the thermal medium in heavy-ion collisions. 
The remainder of the paper is organized as follows: In 
Section II, we review single inclusive hadron fragmenta- 
tion functions in terms of overlapping matrix elements 
between parton field operators and hadronic states. In 
Section III and IV, we introduce the hadronic wavefunc- 
tions in the constituent quark model and reformulate 
the parton fragmentation functions in terms of quark re- 
combination probabilities and shower quark distribution 
functions of the fragmenting partons, which are defined 
as the overlapping matrices of the parton field opera- 
tors and constituent quark states. We also derive the 
QCD evolution equations for the shower quark distribu- 
tion functions and discuss possible sum rules relating tri- 
quark, diquark and single quark distribution functions. 
In Section V, we extend the formalism to parton frag- 
mentation in a thermal medium within the framework of 
field theory at finite temperature and derive various con- 
tributions to the effective fragmentation functions from 
shower-thermal and thermal-thermal quark recombina- 
tion. We summarize the results and discuss future work 
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in Section VI. 



II. SINGLE HADRON FRAGMENTATION FUNCTION 



For hadron production processes, in e + e~ annihilation for example, that involve a large momentum scale, the 
inclusive cross section at leading twist can be factorized into a hard part for parton scattering at short distances and 
a soft part for hadronization at long distances as the parton fragmentation function. Though one can systematically 
calculate the hard part in pQCD due to asymptotic freedom, the parton fragmentation function is nonperturbative 
and currently can only be measured in experiments. However, one can derive the DGLAP evolution equations with 
the momentum scale, which have been tested successfully against the experimental data |30|. In this section, we 
review the parton fragmentation functions as defined in the form of matrix elements of parton operators and the 
corresponding DGLAP evolution equations. 

As in a previous study of medium modification of fragmentation functions with detailed balance [3l| , we consider 
e + e~ annihilation to illustrate the parton operator definition of the quark fragmentation functions. Within the 
collinear factorization approximation, the inclusive differential cross section can be expressed as 



dzh 



-*h 1 e 4 dW ^ ^ ? . 



D*(z h )] , 



(1) 



where q = p a + p^ is the four-momentum of virtual photon and s — q = Q 
leptonic and hadronic tensors are given by L^ v (j> a ,pb) = (l/4)Tr (7^^71/2^) an d 



is the invariant mass of e + e pair. The 



W^(q) = £<0 | J" (0)| X) (X |J" (0)| i)){2^5\ Px q) = f 



= / d*xe-^- x (0|J"(0)J"(aO|0), 



(2) 



where J2x runs over au possible intermediate states and the quark electromagnetic current is J M = eqipq^^ip, 



Here, e q is the fractional charge of the quark in units of an electron charge. The total cross-section at the lowest order 
in pQCD is CTq 9 = N c 4ira 2 eq/3s, where N c = 3 is the number of colors in the fundamental representation of SU(3) 
and a is the electromagnetic coupling constant. The single inclusive fragmentation functions of a current quark and 
antiquark are defined as |32l l33l l34l l35| . 



D h q {z h ) = |/ 



" / dip x ( '» Ph ' n 
2 ./ (2n)^{ Zh p-n 



D h Az h ) 



(27T) 4 



S z h - 



Ph ■ n 
p ■ n 



d 4 xe- ipx Tr 



d 4 xe- ipx Tr 



J2(0 1-0(0)1 S lPh )( Ph , S \i>(x)\ 0) 



j2(o\m\s, Ph : 



7 • n 

2p h ■ n 



(3) 
(4) 



respectively, where ph is the four-momentum of the identified hadron and Zh is its momentum fraction with respect 
to the initial parton momentum. Sum over the color index of the quark field is implicitly implied. The light-like 
vector is defined as — [n + ,n~ ,n±] = [0, 1, OjJ. These fragmentation functions DK-Az^) can be interpreted as the 
multiplicity distribution of hadrons with fractional momenta between z^ and z% + dzh produced in the fragmentation 
of the initial parton quark q or antiquark q. 

In deriving Eqs. (JTJ, © an d |JH from Eq. ©, the intermediate states ^ x \X)(X\ have been replaced by two 
complete subsets s , \S,ph', S')(S'; S,Ph\, corresponding to independent fragmentation of the quark and anti-quark 
created from the e + e~ annihilation. One can further replace \S') by a complete set consisting of quarks and gluons. 
To the lowest order of pQCD, only a single quark or anti-quark state \p') of \S') contributes to the hadronic tensor in 
Eq. J2J). Thus, one can effectively express, 



d 3 P ' 



{2tt) 3 2E p 



\S,Ph;p')(p';Ph,S\ 



(5) 



and the factorized form of the differential cross section in Eq. Q can be represented by a cut diagram as shown in 
Fig. ^ in which the amplitudes of the fragmentation processes are represented by the blob connecting initial quark 
and final hadron states. This is normally presented as cut- vertices |31|. Since the hadronic states \X) are all color 
singlet, the above approximation implies that \S) also carries color which can only be neutralized via exchange of soft 
partons. These soft processes are higher twist and are suppressed by powers of 1/Q 2 . Soft gluon exchanges also lead 
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FIG. 1: diagram for single hadron fragmentation in e e annihilation. 



to eikonal contributions or gauge links that ensure the gauge invariance of the defined quark fragmentation functions. 
They, however, do not appear in the light-cone gauge, n ■ A = 0. 



Similarly, the hadron fragmentation function D^(zh) from a gluon parton is defined as 



D h ( Zh ) = ^ f — 



d i p_ 



S z h 



Ph ■ n 
p ■ n 



d 4 xe~^ x d^ (p) £<0 1 A" (0) | S, Ph ) ( Ph ,S\A"(x)\ 0) 



(0) 



where gluonic color indices are also implicitly summed over. The gluon polarization tensor is defined in terms of the 
polarization vectors e M (k) and has the form 



d» v (p) = ^2 e f _ l (p,X)e^{p,X) 



A=l,2 



Pp,n v + p u n^ 
n ■ p 



(7) 



in the light-cone gauge. 

The scale dependence of the single inclusive hadron fragmentation functions can be calculated by considering the 
radiative corrections to the fragmentation functions at the next-to-leading order in perturbation QCD [33l l34j . It is 
governed by the DGLAP evolution equations [53, IH, , 



Q 2 4^D h M>Q 2 ) = asiQ2) [Uz 



Q 



dQ 2 
d 



2tt 



lqq ( z )D^,Q 2 )+ lqg (z)D h g (^,Q*) 



lgq (z)D*(^,Q 2 )+ lgg (z)D h g (^,Q 2 ) 

z z 



dQ 2 •" v "" " 2;v 
where D^(zh/ z,Q 2 ) is the singlet quark fragmentation function 

D':(z, Q 2 ) ee ]T [D$(z, Q 2 ) + D*(z, Q 2 )} , 



(8) 
(9) 

(10) 



and the splitting functions are |36( 

lqq( Z ) = CF 

Jgg(z) = C F 



l + z 2 3 r/1 , 
(l-z)+ 2 
l + (l-z) 2 



lgq (z) = T F z 2 + (1 - zf 



Jgg(z) = 2C A 



l-z 



{i-zy 



(11) 
(12) 
(13) 
(14) 
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Here, n/ is the number of quark flavors, the SU(iV c ) Casimirs are given by Cf — (N 2 — 1)/2N C , Ca — N c and 
Tp = 1/2. The '+ '-function is defined such that the replacement 

is valid for any function f(z) that is continuous at z = 1. 

The single hadron fragmentation functions as defined in Eqs. J3J), l@J and JfjJ) satisfy the following momentum sum 
rules, 

J2 J dzzD h a (z, Q 2 ) = 1, (a = q, q, g). (16) 

h 

One can also naively define the zeroth moments of the fragmentation functions, 

£ J dzD h a (z,Q 2 ) = N h g , (17) 

(18) 



h 



which are simply the average hadron multiplicities of the parton jets. In principle, the average hadron multiplicities 
are not infrared safe and therefore not well defined in the simple leading log and leading twist approximation. One 
has to go beyond these approximations and take into account the coherence of parton cascade. In practice, one can 
introduce a cut-off to regularize the infrared behavior. 

III. QUARK RECOMBINATION AND PARTON FRAGMENTATION 
A. The Constituent Quark Model 

In hard processes involving a large momentum scale, the produced partons are normally off-shell and subject to 
radiative processes following the hard scattering. The DGLAP evolution equations describe these radiative processes 
for a highly virtual parton in pQCD. Eventually, however, when the parton's virtuality becomes smaller than allowed 
for the applicability of pQCD, further interaction among produced or shower partons and the ensuing hadronization 
can never be described by pQCD. Therefore, one cannot calculate perturbatively the parton fragmentation functions 
with momentum scale Q < Qo- Conventionally, one uses experimental data to parameterize the fragmentation 
functions with Q < Qo (called initial condition) and the DGLAP equations [Eqs. © and 0] predict the evolution of 
the fragmentation functions with the momentum scale. 

In this paper, we will employ a constituent quark model to describe the soft interaction and hadronization of shower 
partons. In this approach, soft interaction between shower partons below the scale Qo will be effectively represented 
by constituent quarks and their interaction. The non-perturbative conversion between shower partons and constituent 
quarks will be described by the constituent quark distribution functions in the fragmenting jet. Further interaction 
among constituent quarks during the hadronization will be given by the hadrons' wave functions in the constituent 
quark model. 

In this constituent quark model, a hadronic state with momentum ph = [p~^ , 0, 0j_] can be expressed as 

\Ph) = J[d 2 k ± ]{dx}<p h (k 1± ,x 1 - 1 k2±,X2)\k 1± ,x 1 - 1 k 2 ±,X2), (19) 

for a meson and 

\ph) = J[d 2 k ± }[dx}if h (k 1±7 x 1 -,k2±,X2;h ± ,x 3 )\k 1± ,x 1 ;k2±,X2;k 3± ,x 3 ) , (20) 

for a baryon. Here, 

(n \ n , 2 , 

kiX II 0(0% ' ( n • = 2 for a meson and 3 for a bar y° n ) > ( 21 ) 
i=l / i=l ^ n ' 

[dx] ee S l-5> 11-^' ^ = ( 22 ) 
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FIG. 2: The cut diagram for single meson fragmentation functions in terms of quark recombination. The blob connecting initial 
quark and final constituent quarks before recombination represents the amplitudes of shower-to-constituent quark conversion 
of the fragmenting quark. Hadronic wavefunction is presented by the shaded blob. 



and <ph is the hadronic wavefunction. In the infinite momentum frame, the constituent quarks can be considered 
approximately on shell. The normalization of the hadronic wavefunction is 

[[d 2 k x }[dx}\<p h (k iU Xi)\ 2 = 1, (23) 

given the normalization of single constituent quark states, 

(kix,X i \k j ±,x j ) = (2n) 3 2x l 5 (2) (kix - k j± )S{xi - x ) . (24) 



B. Double Constituent Quark Distribution Functions 



Substituting the hadronic state in Ea. ll9|) in the constituent quark model into Eq. ©, we can express the meson 
fragmentation as 



D q (zm) = — J j^S [ zm- 



-^j J d 4 xe-' lp - x [d 2 k ± }[dx}[d 2 k' ± }[dx'} 



fM(ki±,xi;-ki±, 1 - xi)ip* M {k' li _,x' 1 ] -k' 1± , 1 - x[) 



Tr 



2p + 



^(0 1^(0)1 S; ki±,Xi;-hi±, 1 - xi)(-k' 1± , 1 - x[; k' 1± , x[; S \ip(x) \ 0) 



M S 



(25) 



Since the constituent quarks can be regarded as a complete set of intermediate states \S), they are equivalent to 
and can replace the complete set of hadronic states \S) in the above expression, J2 S \S)(S\ = ^§ O ne can 

then interpret the above expression in terms of quark recombination as shown in Fig. [2 a quark and anti-quark 
with momentum [x\p M , fcijj and [(1 — x\)p M , — from the fragmentation of the parton jet will combine to form 
a meson with momentum pm — \p M ,0, 0j_]. Since the fragmentation function D^(zm) is proportional to the single 
inclusive cross section for hadron production, it should naturally include the interference between recombination 
of quark and anti-quark pairs with different momenta that sum to the meson's momentum. With the presence of 
the interference contribution, one can never arrive at a probabilistic interpretation as postulated in all the current 
recombination models |24l l25l l26 | . However, if the hadronic wavefunction <fM(ki±, x±; k2±,x<i) in the constituent 
quark model is sharply peaked, one may neglect the interference contributions. Effectively, one can complete the 
integration J[d 2 k' ± ][dx'} in Eq. (|25|l and assume that the final result is proportional to the diagonal term with a 
coefficient Cm- One has then approximately, 



D™{zm) 



W,f — — 



Tr 



d 4 p 
(2tt)4 



2Pm s 



E 

s 



d 2 k^ 



dxi 



(2tt) 3 4a;i(l-o;i) 



\<PM (fax,xi; -kix, 1 - xi)Y 



(0 |V>(0)| S;k 1± ,x 1 ;-k 1 x, 1 - Xi)(-kix, I - x 1 ;k 1± ,x 1 ; S \ip(x)\ 0) 



(26) 
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where Cm is a constant with dimensions of momentum that accounts for the non-diagonal (or interference) contribu- 
tions. Note that 



Pm ~ P + Pm ~ z m 



(27) 



and Zi is the fractional momentum carried by the constituent quark. Defining the diquark recombination probability 

Rm as 



Ru(kix, —) 
zm 



Vm^u, — ; —ki±, 1 -) 

zm zm 



(28) 



we can approximate the meson fragmentation function [Eq. I|2fijl] as 



Df{z M ) « Cm / =tRm(0±, —)F^ (z 1; z m - Zl ) , 
Jo 2 



(29) 



where the double constituent quark (or diquark) distribution function in a quark jet is defined as 



-M 



2z lZ2 



Tr 



2p 



A d 2 k^ 



d 4 P 



/ _J_ X y _ f_M 

2(2^ J {2-kY°\ ZM p+ 



Pm 



-ip-x 



T \ - 



(0 |-0(O)| 5", k u k 2 )(k 2 ,k u SU{x)\0) 



M 



(30) 



Here, pm = k% + k 2 , zm = Zi + z 2 , and A is the cutoff for the intrinsic transverse momentum of the constituent quark 
inside a hadron, as provided by the hadron wavefunction. One can consider A as the scale for hadronization. One can 
also similarly express single inclusive hadron fragmentation functions for antiquarks and gluons in terms of the same 
diquark recombination probability Rm{x\,x 2 ) and double constituent quark distribution functions from antiquark 
and gluon jets, F? iq2 (zi, z 2 ) and F qi ^ 2 (zi, z 2 ), respectively, which are defined as 



Ft 2 ( Zl ,z 2 ) 



F^(zi,z 2 ) 



A ,2 



d k x x_ f d p 



(2n) 



2z 1 z 2 J 2(2?r) 3 

Tr J2 (° F(°) I S, fci , k 2 ) (k 2 , k u S \if>(x) | 0) 



2p 



M 



A .2 



d 2 h± f d 4 p I p M 
S z M x- 



2z x z 2 J 2(2tt) 3 J (2tt) 
d^(p) l^"( )l S, h, k 2 ) (k 2 , ki, S \A"(x)\ 0) . 



(31) 



(32) 



As we will show in the next subsection, the above definitions of diquark distribution functions have almost the same 
form as dihadron quark fragmentation functions [35j j , except the cutoff A for the intrinsic transverse momentum which 
is similar to the cut-off fi± for the intra-jet transverse momentum in dihadron fragmentation functions. 

The quark recombination process in our constituent quark model happens at a scale /x ~ A during the evolution of 
the jet. Therefore, the expression for hadron fragmentation function, Eq. (|29|l . in this model is valid for jets produced 
at any initial scale Q. Since the recombination happens at a fixed momentum scale /i, the recombination probability 
should be independent of the initial scale Q. The scale dependence of the hadron fragmentation functions D^(zm, Q 2 ) 
will be completely determined by the scale dependence of the diquark distribution functions F qiq2 (zi, z 2 ,Q 2 ). Similar 
to the derivation of the evolution equations of the dihadron fragmentation function |35| . one can also obtain the 
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radiative corrections to the double constituent quark distributions from parton fragmentation, 

a s (/j, 2 



F^( Zl ,z 2 ,Q 2 ) = F^( Zl ,z 2 ,^) + 



* U ki J Zl+Z2 z^ {z) ^ \z> z" 



+ 



a s (p- 2 ) 
2ir 

2ir 



F^(z u z 2 ,Q 2 ) = F^(z 1 ,z 2 , f t 2 ) + 



f Q dk 2 ± f 1 dz (z t z 2 



(33) 



2tt 
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Z1+Z2 



lgq {z)F? * (* 
\ z 



z 2 2 



a s (M 2 ) 



dfc 2 



+ 



2tt 


A- 
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« S (M 2 ) 






2tt 




z 2 




f 




2tt 




T 2- 



_L J Z±+Z 2 



( 



Zl ^2 2 
— j — 

z z 



i=l <; 

1-^2 



dz 



2(1-2) 



dz 

The singlet constituent diquark distribution function is defined as 

Ff- *(zi J z 2> /i 2 ) = ^ [Ff *(z 1 ,z 2 , M 2 ) +F| 1 ^(z 1 ,z 2;M 2 )] 



z 2 2 



(34) 



(35) 



The splitting 7-functions are the same as in the evolution equation of single inclusive hadron fragmentation functions 
given in Eas. (|ll[) - I|14|) . The first two terms in Eqs. (|33|l and (|34|l represent diquark production from a single parton 
after the radiative splitting. The third term in Eq. I|33[) and the third and fourth terms in Eq. (|34l) are from 
independent single quark production from each of the partons after the radiative splitting. The indices % = 2, 1 and 
Si = <?2? <7i for i = 1,2 and qi = qi,q 2 , representing the exchange between qi and q 2 in the independent single quark 
production process. These terms from independent single quark production are proportional to the products of two 
single constituent quark distribution functions, 



F q q {z q ) 



ij 


f d±p 








f d^p 




(2vr) 4 




f d A p 




(2tt)* 



-ip-x 



Tr 



7 • n 



-X>|V(0)|£,p g >(p g ,S|#s)|0) 



-ip-x 



Tr 



2rf 



]T<0 |^(0)| S,p,)2-£<p„ S \i>(x)\ 0) 



<5 U - ^ 



(36) 

(37) 
(38) 



whose definitions have the same form as single hadron fragmentation functions in Eqs. ©,® an d ©, replacing the 
single hadron state h with a single constituent quark q. 

In the independent single quark production, there are no virtual corrections. Therefore, the corresponding splitting 
7-functions 



7 99 ( 



C, 



1 + z 2 



% g (z) = 2C A 



1-z ' 

z f — z 



1-z 



z 



+ z(l-z) 



(39) 
(40) 



have no '+'-function and delta functions, unlike "/ qq and "f gg in Eas. (|ll|) - (|14|l . 

Since the hadrons' wavefunctions in the constituent quark model restrict the relative transverse momentum of 
the constituent quarks within a hadron to a finite value, we have an intrinsic transverse momentum cutoff A in the 
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definition of the double constituent quark distributions that are relevant for quark recombination in Eqs. (|30() - (|32|l . 
Similarly, such a restriction should also be applied to the relative transverse momentum between the two quarks 
from the independent single quark production processes in the radiative corrections. This is why there is a cutoff 
A in the relative transverse momentum in the third term in Eq. (|33|) and the third and four terms in Eq. (|34|) for 
the independent production of two constituent quarks that will form a final hadron. This cutoff is determined by 
the hadronic structure in the constituent quark model and is independent of the momentum scale of the initial hard 
parton scattering. Therefore, these radiative corrections do not contribute to the DGLAP evolution equations for the 
diquark distribution functions. Differentiating Eqs. <|33|) and l|34|) with respect to Q 2 , one obtains, 



Q 2 ^F^( Zl ,z 2 ,Q 2 ) 



a s (Q 2 ) 
«,(Q 2 ) 



zi+z 2 

2\ /•! 



dQ 2 



2n 



dz 

dz 
^2 



lgq{z )F^(^r-,Q 2 )+l g9 {z)FI^{^^,Q 2 ) 



(41) 
(42) 



One can show that combining the above evolution equations for the double constituent quark distribution functions 
with the expression of jet fragmentation function in Eq. (|29|l . the DGLAP evolution equations for single inclusive 
hadron fragmentation functions in Eqs. ijBJ and © can be recovered. 

Similar to the single and dihadron fragmentation functions, one can also derive sum rules for the single and diquark 
distribution functions in a parton jet. One can define the mean constituent quark multiplicities from the single quark 
distributions as 



dzF2(z) = Nl(a = q,q,g). 



Similarly, single quark distribution functions also obey the momentum sum rule, 

/ 1 dz^zF ( ?(z) = l. 
Jo 



(43) 



(44) 



For diquark distribution functions F^ 1,q2 {z\ 1 z 2 , Q 2 ), the multiplicity sum rule leads to the second cumulant moments 
of the multiplicity distribution, 



J d Zl dz 2 F^{ Zl ,z 2 ,Q 2 ) 



2^ — 



(45) 



where S qiq2 = 1 if q\ and q 2 are identical quarks and 6 qiq2 = if q\ and q 2 arc different. It illustrates the importance 
of quark correlations in the diquark distribution functions. Any initial conditions for the diquark distributions should 
contain the correlation between two quarks within a jet. For two identical quark distributions, the above equation 
also illustrates the information on the multiplicity fluctuation contained in the diquark distributions. Because of the 
correlation and fluctuation contained in the diquark distributions, one cannot find any rigorous sum rules relating 
single and diquark distributions. Following the work on dihadron fragmentation functions 35], an approximate ansatz 
that contains the minimum correlation and fluctuation is 



J2 [ dz 2 F^( Zl ,z 2l Q 2 ) 

92 J 



Nl 1 



(46) 



Using the momentum sum rule Eq. Ij44fl for single quark distributions, the above ansatz will lead to the momentum 
sum rule for double constituent quark distributions, 



f d Zl dz 2 ^( Zl + z 2 )F^( Zl ,z 2 ,Q 2 ) 

9192 



Nj 1 (Net 2 ~ 5 qi ^ q2 ) .— gi "5-^92 s 
-91 -irr-1-2 V iy a ~r Jv a 1 ' 



(47) 



Whether such ansatz is a good approximation remains to be explored. 



C. Dihadron Fragmentation Functions 



The definitions of diquark distribution functions of a jet are very similar to that of dihadron fragmentation functions, 
except that the constituent quarks in the distributions are restricted by the hadron wavefunction in the phase space 
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for quark recombination during hadronization. According to Ref. 35], the dihadron fragmentation functions of a 
quark or gluon parton are defined as, 



D^ h *(z hl ,z h2 ) 



D^(z hl ,z h2 ) 



D^(z hl ,z h2 ) 



d 2 Ph 1± f d 4 p 



2z hl z h2 J 2(2^)3 J (2iry 



S ( zh - 



pi 



Tr 



7T~f }_^( Q \^( Q )\ s ,Phi,Ph 2 ) (Ph 2 , Phi , S | ip(x) | 0) 
2 Ph s 



d 2 Ph 1± f d 4 p 



2z hl z h2 J 2(2^) 3 J (2 



-ip-x 



Tr 



5>|#))|£,p &1J p fca ) ^ + (Ph 2 , P hl ,S\i>(x)\0) 

S Zp h 



! d 2 p hl± f dp i _p h 
r° \ z h , 



2z ftl z ha J 2(2tt)3 J (2tt) 
V (?) E (° I ^ (°) I 5 ' Phi » P^ ) > Phi . S 1^ (*) I °) > 



(48) 



(49) 



(50) 



where ph = Ph x + Ph 2 an< i z ft = z h t + ^h 2 - Note that the transverse momentum, Ph t ± and Ph 2 ±, are defined as 
perpendicular to the summed total momentum p^. The relative transverse momentum used in Ref. |35| is then 
q± = 2p hl± = -2p h2± . 

The radiative corrections to the dihadron fragmentation functions are similar to that of diquark distribution func- 
tions, except that the transverse momentum between two individual hadrons is not limited as in the case for diquark 
distributions due to quark recombination during hadronization. Therefore, the cutoff A in Eq. Ij33|l and (134ft should 
be replaced by Q 2 . The DGLAP evolution equations for dihadron fragmentation functions then become |35| . 



D h g ^(z hl ,z h2 ,Q' 2 ) 



dQ 



2 3 



2ir 

2 

+E 

1 = 1 JZ K 



Jz hl +z h2 z \ Z Z 



a s (Q 2 
2tt 

2 



Zh-\ +Zh 



dz 

72 



EE 



dz 



., z(l - z) 

1=1 q " z hi v ' 



(^(^,Q 2 )d^(^-,q 



dz 



l gg {z)D^(^f, Z -f,Q 2 ) 



dz 



2(1 -Z) 



lag 



(51) 



(52) 



The additional terms in comparison to the DGLAP evolution equations for single hadron fragmentation functions that 
are proportional to the convolution of two single hadron fragmentation functions arise from independent fragmentation 
of both the leading parton and radiated gluon. They can be attributed to the dependence of the relative transverse 
momentum between two hadrons from the independent fragmentation on the momentum scale. This is the main 
difference between the evolution of the above dihadron fragmentation functions and the diquark distributions in the 
quark recombination model. 
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IV. BARYON FRAGMENTATION FUNCTIONS IN QUARK RECOMBINATION MODEL 
A. Triple Constituent Quark Distribution Function 



Since a baryon has three constituent quarks, the fragmentation functions for a parton into baryons naturally involve 
triple quark distribution functions in the quark recombination model. Similar to the meson fragmentation functions, 
we can also express the baryon fragmentation functions in terms of baryon wave functions and overlapping matrix 
elements of parton fields and constituent quark states, as illustrated in Fig. |3 



Df(z B ) = 



d 4 p 



5 I zb- 



^\ J d 4 xe~ tp - x [d 2 k ± }[dx}[d 2 k' ± }[dx'} 



<PB(ki±,xv, k 2 ±,x 2 ]k 3 _ L: x 3 )ip* B (k[ ± , x[; k' 2± , x' 2 ; k' 3± ,x 3 ) 



Tr 



y^(0 |^(0)| S; kx±, n; k 2± ,x 2 ; k 3± , x 3 ) {k' 3± ,x 3 ; k' 2± , x 2 ; k' 1± , x[; S \tjj(x) | 0} 



Or — 



d 4 p 

WT 4 



5 z B 



-Eb\ f -ip-x f d 2 k 1± dxi f d 2 k 2 ± dx 2 1 
P+J J C ' J (2tt)3 AxJ (2tt)3 4x 2 x 3 l 



\<fib(ki±,xi;k 2± ,x 2 ; k 3J _, x 3 )Y 



Tr 



2pj 



y^(0 |^(0)| S; kix, xi\ k 2± ,x 2 ; k 3± , x 3 ) (k 3± ,x 3 ; k 2± , x 2 ; k x ±, x\\ S \tp(x) | 0) 



(53) 



where [d 2 fcj_] and [dx] are given by Eqs. 12 1H and (12211 for baryons, x 3 = 1 — x\ — x 2 , k 3 ± = — fcu_ — k 2 j_, and C'b is a 
constant with the dimension of momentum. Notice again that 



X\ t - ~r !I7 ? X 2 

Pb P Pb Zb 



z 2 



(54) 

Pb P Pb Zb 

As in the quark recombination model for meson fragmentation functions, we also define the recombination probability 
for a baryon as 



-Rs(fci_L, — ; fc 2 j_, — 
zb zb 

The baryon fragmentation function from Eq. 153f) can then be cast in the form 



, n Zl Z 2 Zl+Z 2 2 

\<Pb(ki±, — ; fc 2J _, — ; -feu. - fe 2 ±, 1 ) 

zb zb zb 



d . . „ r ZB d Zl r- 

D*(zb)*CbJ q ~yJ o 



2 z B z B 



(55) 



(56) 



This relation is very similar to that for meson fragmentation functions in Eq. I|29|) . The triple constituent quark 
distribution function is defined as 



A d 2 k 1± _ f A d 2 k 2 ^ 



d 4 p 

2 Zl z 2 z 3 J 2(2tt) 3 J 2(2tt) 3 J (2tt) 4 

7^ 



Tr 



; f ^(0 \i>{0)\S, k 1 ,k2,k 3 ){k 3 ,k 2 ,k u S\i>{x)\0) 
1 Pb ~ s 



(57) 



where, ps = ki + k 2 + k 3l zb = Zi + z 2 + z 3 and A is the cutoff for the intrinsic transverse momentum of the constituent 
quarks in a baryon. Similarly, the triple constituent quark distribution function from the antiquark and gluon parton 
are defined as 

rA d 2 k 1± f A d 2 k 2 



{Z\,Z 2 ,Z 3 ) 



d 4 p 

2 Zl z 2 z 3 J 2(2tt) 3 J 2(2tt) 3 J (2tt) 4 ' 



5 zb- 



Pb 
p+ 



Tr 



^(0 |^(0)| S, k x> k 2 , fta)-hjr(A;a, fc 2 , k u S \^(x)\ 0) 

a ^PB 



Df^( Zl ,z 2 ,z 3 ) = 



A >2 



d z k 1A 



2z lZ2 z 3 J 2(2tt) 3 J 2(2tt) 3 7 (2tt) 4 



A d 2 k 



2± 



d 4 p 



3 ( zb 



Pb 



d i xe~ lv - x 



d^(p) J> |^(0)| S, k u k 2 , k 3 )(k 3 , k 2 ,h,S \A»(x)\ 0) 



(58) 



(59) 
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Similarly as for the diquark distribution functions, one can also obtain the radiative corrections to the triple 
constituent quark distribution functions, 



F^( Zl ,z 2 , z 3 , Q 2 ) = F^( Zl ,z 2 , z 3 ,/i 2 ) 

2tt 7^2 kj_ J Z1+Z2+Z3 z 3 q \z z z ) 

^/;f/L +! ,l-W»(f.?.?.^) 

^eCHC"^*"* cm*- fees) 

( Z1 ,Z2, za, Q 2 ) = ( Zl ,z 2 , z 3 , M 2 ) 

^/;f/L + „i-<^-(?.?.?- s ) 



where (i', i") — (2, 3), (1, 3), (1, 2) as i — 1, 2, 3 and the singlet triple constituent quark distribution is given by 



Fr 293 (zi, z 2 ,z 3 , M 2 ) = ]T [F^(z u z 2 , z 3 , (j?) + F^(z u z 2 ,z 3 , M 2 )] 



(62) 
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The DGLAP evolution equations for triple constituent quark distribution functions can be obtained from the above, 



Q 2 ^ 2 -Df™°(z 1 ,z 2 ,z 3 ,Q 2 ) 



as(Q 2 
2tt 



dz 
Z1+Z2+Z3 



Z\ Z2 Z 3 



7 99 (^ 19293 (-,-,-,Q^) 



Z z z 



z\ z 2 z 3 



z z z 



a s (Q 2 ) 



dQ 



2tt 



Z1+Z2+Z3 

Si 

9 



7«,(*)tff lfc *(-,-,^,Q 2 ) 

z z z 



(63) 



(64) 



Again, the radiative corrections from the independent fragmentation of the leading and radiated partons do not con- 
tribute to the DGLAP evolution because of the restriction on the relative transverse momentum between constituent 
quarks by the baryon's wavefunction during quark recombination. As in the case for diquark distribution functions 
and meson fragmentation functions, the above DGLAP evolution equations for triple quark distributions will give 
rise to the DGLAP equations for baryon fragmentation functions within the quark recombination model as given by 
Eq. JU. 

One can also obtain the multiplicity sum rule for the triple quark distribution functions of a parton jet, 



dz x dz 2 dz 3 F^ {zi,z 2 , 23, Q ) = Nl 1 {N q a 2 - 5 qiq2 ){N q a 3 - 5 q2Q3 - 6 qiq3 ) 



(65) 



Similar to the ansatz of the sum rule relating single and diquark distributions, one can also have the following ansatz 
for a sum rule relating diquark and triple quark distributions, 



£ J dzsFr 2q3 (zi,z 2 ,z 3 ,Q 2 )^ ^^ AfiPh 5q : q3 w fr 2 (^ 2 ,Q 2 ). 



m i (N q a 2 -S qiq2 ) 



(66) 



B. Trihadron Fragmentation Function 

It is straightforward to generalize the result for the triple quark distribution functions to the trihadron fragmentation 
functions. The operator expressions of the trihadron fragmentation functions from quark, antiquark or a gluon jet 
are defined as 



D 



h 1 h 2 h 3 



(Zfn , Zh 2 , Zh 3 ) — 



D^ h2h3 (z hll z h2 ,z h3 ) 



J^hih2h3 



{ z h± 1 , Zfr 3 ) 



d 2 Ph 1± f d 2 p h . 



d 4 p 



S ( zh- 



Pu 



2z hl z h2 z h3 J 2(2tt)3 7 2(2^)3 ,/ (2tt) 4 V 

+ _ 

Ph 2 ,Ph 3 )(Ph 3 ,Ph2,Phi,S \ip(x)\ 0) 



Tr 



z Ph s 



d 2 Ph 1± f d 2 p fl2 



d A p 

2z hl z h2 z h3 J 2(2tt) 3 J 2(2^) 3 J (2tt) 4 



Tr 



2^(o 1 ^(0) 1 s, p hl , p h2 , Ph 3 ) — - (ph 3 , Ph 2 , Ph x , s 1 t/j(x) 1 0) 

c Z Ph 



d 2 Ph 1± f d 2 p h . 



d 4 p 



* ( *h - 4 



d i xe- %px 



2z hl z h2 z h3 J 2(2n) 3 J 2(2^)3 J (2tt) 



(67) 



(68) 



(69) 



where p h = p hl + p; l2 + ph 3 and z h = z hl + z fl2 + z fl3 . 

One can show that in the leading twist and collinear factorization approximations, the semi-inclusive cross section 
°e+e-— yh!h2h 3 can be expressed in terms of the trihadron fragmentation functions as, 



da, 



e+e — >h\h2hi 

dz hl dz h2 dz h3 



h * h ° (z hl , z h2 ,z h3 ) + h ^ (z hl , z h2 , z h3 ) 



(70) 
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Again, the radiative corrections for the trihadron fragmentation functions are similar to that of triple quark dis- 
tribution functions in Eas. i|(jU|) and (|61|) . except that the cutoff in the intrinsic transverse momentum A should be 
replaced by Q which is the only limit of the relative transverse momenta between hadrons. The DGLAP equations 
for the scale evolution for the trihadron fragmentation functions are then, 



qq 2 D q 1 2 3 ( Z hn z h 2 , z h 3 i Q )— ^ 



Jz hl +z h2 +Zh 3 z - 



E 



1=1 JX *i 

3 „l-z 



~ !: 1 '' " dz , % ,, (z h \ 



E 



z(l - zf 
dz 



a 



Zh { , Zh,„ 



z 2 (l- z) 



dz < ^h lh2 h 3 f z hi z h 2 z h 



1 -z' 1 -z 



a 



z z z 
1 

Zhi + z h 2 + z h 3 



dz 



z z z 



1— Zh., -Zh.,, 



EE 

i—l q 

3 

i=l q Jz h i ,+ Z h^ 



dz 



Zh i 



z{l- zf 
dz 



1m 



(z)D^(^,Q*)d^»(^,^ 



z 



dz 



z h 1 + z h 2 + z h 3 



3 

E 



l-Zh., -Zh.„ 



" z r„\r>Vv f Zh *' Zh -" nA n h i ( Zhi n'< 



dz 

z{l~zf 



Zh z , z h .„ 

T~z'Y~z' 



(71) 



(72) 



V. QUARK RECOMBINATION AND JET FRAGMENTATION IN A THERMAL MEDIUM 

So far we have reformulated jet fragmentation in vacuum in the quark recombination model in which we related the 
parton fragmentation functions to the quark recombination probability as determined by the hadron wavefunction 
in a constituent quark model and multi-quark distribution functions of a parton jet before hadronization. In the 
following sections we will extend the formula to the case of jet fragmentation in a thermal medium which is relevant 
to jet fragmentation in the environment of high-energy heavy-ion collisions. Such problems have been considered 
before [31|. But the attention has been focused on parton emission and absorption by the propagating parton jet in 
a thermal medium before hadronization. It has been assumed that the thermal medium is in a deconfined phase so 
that partons from the jet will eventually hadronize together with the medium. 

In this paper, however, we focus on the physical process during the hadronization of the parton jet in a thermal 
medium and the modification of the parton fragmentation functions with respective to that in the vacuum. We assume 
that the effective degrees of freedom can be described by constituent quarks just before and during hadronization. 
Therefore, we consider constituent quarks not only as the effective states in the process of jet fragmentation but 
also the effective constituent of the medium just before hadronization. In this framework, one naturally encounters 
recombination between thermal constituent quarks and shower quarks from parton fragmentation, in addition to 
recombination of shower quarks as in the vacuum. They both contribute to to hadron production associated with an 
energetic parton jets in a thermal medium. 

As in the previous study [3 1| . one can describe the fragmentation of a parton jet in medium simply by replacing the 
vacuum expectation in the S matrix of the processes or the operator definition of the parton fragmentation functions 
by their thermal expectation values, (O|0|O) — > ((C)), 



(73) 
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where, H is the Hamiltonian operator of the system and 1//3 = T is the temperature. Therefore, the single hadron 
fragmentation functions at finite temperature for a quark, antiquark and gluon jet can be defined as 



D h q (z h ,p + ) 
D h g (z h ,p+) 



d 4 p 

(2tt)* 

d 4 p 

(2n) 4 

d 4 p 

(2^) 4 



pi 



S Zh 



d 4 xe- ip - x Tr 



d 4 xe- ip - x Tv 



z Ph s 



^mo)\s,Ph)f+(ph,s\fi>(x))) 

S ZP h 

d 4 xe- ip - x d^{p)J2(( Afl (0)\S,p h )(p h ,S\A^x))}, 



(74) 
(75) 
(76) 



where ph and p are the four-momentum of the hadron and the initial parton, respectively. In the above definition of 
the parton fragmentation function in medium, we have explicitly kept the dependence on the initial parton energy 
p + . Such dependence arises at finite temperature from the the thermal average which also introduces dependence on 
the temperature T. An alternative definition of the fragmentation functions |3ll |. 



D h M,p + ) 



D h M>p + ) = - 



pi 



d 4 k 



%>+ J (2tt)< 



<5(fc+ -p H 



d A xe- tk - x 



Tr 



- &/t^(* + -p + )/**° 



d 4 k 



Ap+ J (2tt) 4 
d 4 k 



-ik-x 



Tr 



j +J £mo)\s,Ph){p h ,s\i>(x))) 

s 

J2((m\S,p h h + (p h ,S\i;(x))) 



Ph 

2p+ 



(2n) 4 



6(k+-p + ) / d 4 xe- ik - x J2{{F+»(0)\S,p h ){p h ,S\F+(x))), 



(77) 
(78) 
(79) 



explicitly takes into account the dependence on the absolute initial parton energy and the medium temperature. The 
two expressions can be related via the identity, 



p pi ' 



(80) 



At zero temperature T — 0, the above fragmentation functions will be reduced to the parton fragmentation functions 
in vacuum as defined in Eqs.10, l@J and ©. At leading twist, they depend only on the scaling variable, Zh — p~^/p + , 
i.e., the ratio of the hadron and parton energies. 



A. Meson Production from Thermal Quark Recombination 
From Eq. 174J1 , the meson fragmentation function at finite temperature from the quark parton can be expressed as 



Tr 



4 xe~ ip ' x 

1 J (2tt) 4 " V~ JW P + J ' 
[d 2 k±] [dx] [ [d 2 k' ± ] [dx'\ ip M {fa± , X\ ; k 2± , x 2 )ip* M {k[ ± , x' t ; k' 2± , x' 2 ) 

7 H 



2Pm 



2«V(0) |5; k lx , xi ; k 2± , x 2 ) (k' 2± , x' 2 ; k[ x , x\ ; S\ip(x)) ) 



(81) 



To evaluate the thermal average in Eq. QSlfl . it is convenient to use the finite volume quantization (FVQ) in a cubic 
box with length L. Employing FVQ, the momentum integral, ^-function and continuous states can be replaced by 
following summation, Kronecker functions and discrete states, respectively. 

27r\ 3 ^ ^ ^ ,,2Tm x 2nn y 2irn z 

T 2^ 2^ 2^ A— r~ »— r - >—r~ )> ( 82 ) 



dq x \ dq y / dq z f(q x ,q y ,q z ) 

) J — go J — oc 

(2n) 3 S(q x - q' x )5(q y - q' y )S(q z - q' z ) 

1 



L ) ^ ^ ^ J v L ' L 

n x — — oo n v — — oo n z — — oo 



^/2E~p 



I?) 



(83) 
(84) 
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Using Eq. (|73|l , we can express the thermal average of the matrix element in Eq. (|81|l as 

((ip(0)\S; k 1± ,x 1 ; k 2± ,x 2 )(k' 2± ,x 2 ; k' 1± ,x[; S\ip(x))) 

E ({{qi,TH)}\e- fi6 il>(0)\S\ h;h) (k 2 ; k[; S$(x)\{(q u m)}) 
_ {M 

E ({(Qi,ni)}\e-^\{(q hni )}) 



(85) 



where n$ = 0, 1 (for i = 1,2 - ■ ■ , oo) since we have assumed the thermal medium to consist of constituent quarks that 
obey Fermi statistics. The denominator in Eq.JHU can be reduced to 

— Q V E °° 

{Hi} {«,} i=l 

Taking into account that the intermediate states of constituent quarks \kx), \k 2 ), \k[), \ k' 2 ) can contract both with the 
field operators ip and ip and the constituent quark states !{(<&, Tii)}) in the thermal medium, we get 

({ip(0)\S; kix, xi; k 2± , x 2 ) (k' 2± , x' 2 ; k[ ± , x' t ; S\ip(x))) 

- ffl(' + . E-^"* 

\*=i / W 

[((<?i,«i); • • • ; (q oo ,n x )\^(0)\S;k 1 ;k 2 )(k 2 ;k[-,S\'tp(x)\(q o<: ,n oo ); ■■■ ;(<?i,th)) 

+£((?i> n i); " ' ' • • • l"0(O)|5;fei; A; 2 >(fc2; fcj.; 5*1^(^)1 • • • ; {qj,nj); • • • ; (<7i,th)) 
W 1 1 1 1 

+£((91, " ' ;(*'"-«);••• 1^(0)1^; kx;k 2 )(k' 2 ; k{; S\i/j(x)\ ■■■ ; {qj,rij);- ■ ■ ; («i,ni)) 
W 1 1 1 1 

+ £ ('*' ;(«) "/);••• l^(0)|S; kx;k 2 )(k' 2 ;k'x;S\ip(x)\--- ; (g TO , n m ); • • • ; (tfr, n^); • • • ) . (87) 
- , i i i i J 

' i i i i 

For simplicity in this paper, the states \k±), \k 2 ), \k[), \k 2 ) without the connecting lines are meant to be only contracted 
with the field operators ip and ip. The states \(qi, ni)) with the connecting lines are meant to contract only with the 
final states \kx, 2 ) and \k[ 2 ) but without contraction with the field operators ip and if). The states \(qi,Tii)) without 
the connecting lines can still contract with the field operators ip and ip but will not contract with the final states \ki t2 ) 
and \k[ 2 ). 

Let us for a moment neglect the contraction between the thermal states \(qi,rii)} and the field operators ip and ip. 
Using 

((<7i,«i);- • • ; (<?i-i,"i-i); (qi+i,n i+ x)\ ; (q j+ x,n j+ x)\ (q^x^j-x); • • • ; (gi,7ii)) = 5y , (88) 

one can obtain from Eq. (|87|l 

((ip(0)\S; ki±,x 1 ;k2±,x 2 ){k! i± ,x' 2 ; k' 1± ,x' 1 ;S\ip(x))) 

n(l + e-^) [£e ^ \0\m\S-,kx-,k 2 )(k 2 ;k'x-,S\iJ(x)\0) 
i=l / {m} 

, -/3E"* E ' ~ ~_ 

3 {",} 1 1 1 1 

,/ 
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+E e "^ E' e 4=1 (fe>^)lV-(0)|5; kx;k 2 )(k 2 ;k'x; SVP{x)\{q 3 ,n 3 )) 
3 {m\ ' 1 1 1 



+ ^ e -/3(E 3 +E m )^ e i=1 (( % ,n J );(< Zm ,n m )|^(0)|5;fc 1 ;fc 2 )(^;fc; ; ^|^(x)|( (Zm ,n m );fe,n,))J ,(89) 
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FIG. 4: The cut diagram for contributions from (a) shower-thermal and (b) thermal-thermal quark recombination to the single 
meson fragmentation function. 



where E , S indicates that i ^ j and i ^ j,m in their summation, respectively. One can generalize the above results 
to include the contraction between thermal states and the quark field operators. In this general case, the vacuum 
expectation will be replaced by a thermal average without contraction between thermal states and the final constituent 
quark states. Using Ea. H86|) and replacements l|82l) . H84[) . we obtain the thermal average value of the matrix elements 
in the meson fragmentation function as 

«V(0)|S; fc ix , Xl ; k 2± ,x 2 )(k' 2± , x' 2 ; k' 1± , x' i; S\^(x))) = fc «V(0)|5; h; k 2 )(k' 2 ; k[; S$(x))) k 
d 3 qi 1 



{2n) 3 2E qi e P E "i + 1 
d 3 qi d 3 q 2 



fc ««iM0)|S; fci; k 2 ){k' 2 ; k[; S\i>(x)\qi)) k + fc ((gi|^(0)|S; kv, k 2 ) (k' 2 ; k[; S\^(x)\ qi )) k 



(2ir) 3 2E qi (2ir) 3 2E q2 eP E *i + 1 e 0E ^ + 1 



fc«<&; qi\1>(0)\S\ kv,k 2 )(k' 2 ; fci; S\^(x)\ qi ; q 2 )) k , (90) 



where the notation &((•••))& represents thermal averaging without the contraction between thermal states and the 
final constituent quark states. The form {e@ E + is the Fermi-Dirac thermal distribution function in the co-moving 
frame of a fluid that maintains only local equilibrium. In the reaction frame in which a parton jet propagates, it 
should be replaced by 

/(«) = ^ri> ( 91 ) 

for the thermal quark distribution in a fluid element with flow four-velocity u in the thermal medium. The momentum 
of a thermal constituent quark is decomposed to q = (q + , 0, q±) with respect to the direction of the propagating parton 
jet. In the following, we also sometimes express the longitudinal component of the thermal momentum as a fraction 
of the initial parton momentum, q + — z q p + . In this case, the value of z q is not bound. With this, we can finally 
express the thermal averaged matrix element in the fragmentation function as 

M0)\S] fcix.xi; k 2X ,x 2 )(k' 2± , 4; fei x ,si; S®(x))) = fe «|V(0)|S; fc x ; k 2 )(k' 2 ; k[; 

' h,1> %f(qi±,qt)\k((qi\m\S; fci; k 2 )(k' 2 ; k[; S\^(x)\ qi )) k + k ({qi\^(0)\S; fc i; k 2 ) (k' 2 ; k[; S\^(x)\ qi )) k 



2(2^)3 q l 
d 2 q 1± dq+ f d 2 q 2± dq. 



2(2^)3 q + J 2(2^)3 q } 



-f(qi±,qt)f(q2±,qt)k((q2; qi\ip(0)\S;kf,k 2 )(k 2 ; k[;S\ip(x)\qi;q 2 ))k ■ (92) 



The above equation has three distinct contributions. One can immediately identify the first term as the contribution 
from recombination of two shower quarks from parton fragmentation, 

b^ ss \z M ) « Cm / ^Rm(0±, —)F^ ( Zl ,z M - zi) , (93) 

l Zm 



(I 



This contribution is normally referred [24J to as "shower-shower" quark recombination. It has exactly the same form 
as the parton fragmentation functions in vacuum [Ea. l|29|) ] in the framework of quark recombination, except that 
the diquark distribution functions F^ 1 ^ 2 (zi, z 2 ) are now medium modified. Their definitions are similar to that in 
vacuum in Eq. (|30|l but the vacuum expectation is replaced by a thermal average. These modified diquark distribution 
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functions should in principle contain effects of multiple scattering, induced gluon radiation and absorption, in the 
same way as the modification of hadron fragmentation functions in a thermal medium |3lj . We have assumed that 
the meson wavefunction and therefore the recombination probability during the hadronization of the parton jets in a 
quark-gluon plasma are the same as in vacuum. 

In the second term, one of the constituent quarks is contracted with the thermal quark states while the other with 
the quark field operators. This corresponds to processes in which a quark from the thermal medium combines with 
another quark from the parton fragmentation (or shower) to form the final meson, as illustrated in Fig.QJa). Using 
the convention of previous studies |24| , these processes are called "shower-thermal" quark recombination. Finally, the 
third term where both of the final constituent quark states contract with the thermal quark states corresponds to the 
formation of the final meson from two thermal quarks in the quark recombination model, as shown in Fig.^Jb). This 
is referred to as "thermal-thermal" quark recombination. 

Using Eqs. (|28|) and l|81l) and the contraction between two single-particle states, 



(qi\kj) = (qix,x qi \kjx,Xj) = (2ir) 3 2x i 6 ( - 2) {qi± - k jA _) S (x qi - xj) , (94) 
i i 

we can express the contribution from "shower-thermal" quark recombination to the meson fragmentation function as 

JO Z M J 2{2iry Z M {l-Zq/ZM) 

f qi (q±,ZqP + )Fq 2 ( Z M - Z q ) + fq 2 {qj_,Z q p + )F^(z M - Zg) , (95) 

where x qi — z q JzM is the thermal quark's momentum as a fraction of the produced meson's momentum, F^(z) and 
(z) are single constituent quark or anti-quark distributions of the fragmenting parton jet in a thermal medium 
defined similarly as in Eqs. (|36f) and l|37|) . except that the vacuum expectation values are replaced again by thermal 
averaged expectation. They should be different from the corresponding quark distributions in vacuum because of 
multiple scattering, induced gluon bremsstrahlung and parton absorption. 

The contribution due to "thermal-thermal" quark recombination to the meson fragmentation function in Eq. I|92|) 
comes from recombination of two thermal constituent quarks. This contribution would never be associated with the 
parton jet had we not used the parton's momentum as a reference to calculate the distribution of produced hadrons 
from such thermal quark recombination. Again, using Eq. I|94(l . we can express this "thermal-thermal" contribution 
as 



f ) M(TT) { 



ZM,P 



x 4 



d 2 q 1 ±dz qi r d 2 q 2A _dz, 
2(2tt) 3 J 2(2vr) 3 

/ Zgl Z q 2 \ 

tpM[qi±, — ;<72_l, 



-fqi (9i-L> z qi p + )f q2 (q 2 ±, z q2 p + )S q {z M ) 



ZM 



ZM 



(27t) 3 5 (2) (pm± - qi± - q2±)S{z M - z qi - z q2 ) 



(96) 



where 



S° q (z M ) 



-M 



d 4 p 



S z M - 



p+ 



d 4 xe- ipx Ti 



52k{(M0)\S){S\l>(x)\0)) k 



ZM 1_ 

2 (2tt) 



d 2 p± 



1 



1 



2(2^)3 z M 



d pm a 



(97) 



is just a factor associated with the phase-space integration. Since the fragmentation function represents the total 
yield of particle production, the thermal contribution naturally contains a volume factor, which leads to the 6- 
function squared in Eq. (|96|l . To factor out this volume factor, we again employ the finite volume quantization with 
V = L? in a cubic box. The momentum integral and (5-functions can be expressed by using replacements Eqs. I|82(l 
and H83f) . Then we can rewrite Eq. I|96|l as 



bf TT \ ZM ,p+) = v P + 



d 2 PM± 
(2tt) 3 



dz„ 



2z M J {2nf 



f qi {q±,z q p + )f q2 (p M ± - q±, {z M ~ z q )p + )R M (q±, — ) • 

ZM 



(98) 



Since hadron production from the thermal quark recombination is not correlated with the parton jet and its 
fragmentation, the above expression is a little unnatural. One should be able to rewrite it in a form that has no 
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dependence on the parton jet. Considering, 



dN M(TT) dN M(TT) 



dzA 



-P 



one can obtain the invariant hadron spectrum from thermal quark recombination, 

d 2 q A 



(2-) 3 f^ =V dx q 

dviid z VM^ Jo 



2(2tt) 



rfgi (q±,XqPtf)fq a (PM± - qx, (1 - x q )pl I )R M (q±,x q ) . 



(99) 



(100) 



which coincides with results from other recombination models |25 . 26] . In this expression, the hadron spectra from 
thermal quark recombination are not correlated and therefore do not depend on the parton jet fragmentation. Even 
though one can sum over the three contributions and obtain the effective meson fragmentation in a thermal medium, 

(101) 



Df (z M ,p + ) = D^ ss \zm) + D^ ST Hz M ,p+) + D^ TT \z M ,p+) , 



the last term from thermal quark recombination is not correlated with the initial parton jet and therefore should 
not be considered as part of the medium modified jet fragmentation function. The contributions that are correlated 
with the initial parton jets are from "shower-shower" and "shower-thermal" quark recombination. For a thermalized 
medium, the thermal quark distribution follows a Fermi-Dirac form that is determined by the local temperature 
and flow velocity. The contribution from "shower-thermal" recombination is most important for hadron spectra in 
the intermediate transverse momentum region. This contribution, however, will be negligible relative to the "shower- 
shower" recombination which dominates hadron spectra at large transverse momentum because that the power-law-like 
spectra of initially produced partons will win over the exponential-like distribution of thermal quarks. 



B. Baryon Production from Thermal Quark Recombination 



We can similarly generalize the quark recombination model for baryon fragmentation functions to the case in a 
thermal medium, 



bi 3 (z B , P +) 



[d 2 k ± ][dx] 



Tr 



Xly/X2y/X3 



[d k' ][dx'] 

f 7= r-r ^B (ki_L, x^, k 2 ±, x 2 ] k 3 ±, x 3 )tp* B (k' 1± , x[; k' 2± , x' 2 ; k' 3± , x' 3 ) 



X -| a / X ^ \/ X 



TTT ^((^(0)15"; fci ± , xi; /c 2 j_, a; 2 ; fax, x 3 ){k' 3± , x 3 ; k' 2± , x' % \ k\±,x' x \S\^{x))) 



(102) 



As in the case of meson fragmentation functions, the thermal average of the matrix element in the above baryon 
fragmentation function can be expressed as 

((ip(0)\S; ki±,Xf, k 2± ,x 2 ; k 3± ,x 3 )(k' 3± ,x 3 ; k' 2± ,x 2 ; k' ix ,x' x ) S\ip(x))) 
= fc «V(0)|§; fci; fo; k 3 )(k' 3 ; k' 2 ; k[; S$(x))) k 

^vTTI -^r/(gl-L)gi")[ k((qi\ip{0)\S; ki;k 2 ;k 3 }(k 3 ;k 2 ;k[;S\iJj(x)\qi)) k 
+ fc ((gi|^(0)|S; h; k 2l k 3 )(k' 3 ; k' 2 ; k[; S^(x)\ qi )) k + fc <( 9 i|^(0)|S; h; k 2 ; k 3 ) (k' 3 ; k' 2: k[; S$(x)\ qi )) k 

I I I I I I I I 

<HhL 1 ,h ' J2L (/ ' i2+ f(qi±, qt)f(Q2±, 91 W0)|S; h; k 2 ; fc 3 ><% fcj ; S|?(af)|fli; 



2(2vr) 3 q + J 2(2tt) 3 q\ 



+ k{{q2]qi\i>^)\S]k 1 ]k 2 ]k 3 ){k' 3 \k' 2 ]k' 1 -,S\^{x)\qi,q 2 )) k + k {{q 2 ]qi\^{Q)\S]k 1 ]k 2 ]k 3 ){k' 3 \k l 2 ]k' 1 - 1 S\^{x)\qi,q 2 )) k 



d 2 qi±_ dqf f d 2 q 21 _ dq£ f d 2 q 31 _ dq^ 



2(2tt)3 q+ J 2(2tt)3 q+ J 2(2tt)3 q+ 

f(qi±,qt)f{q2±,q^)f(q3±,qt) k((q3;q2;qiH(0)\S; k^k^k^^k'^k'^S^^q^q^q^k . 



(103) 
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Here, one can similarly identify the first term and the last term as the contributions from "shower-shower-shower" 
and "thermal-thermal-thermal" quark recombination, respectively. The others terms are the contributions from the 
"shower-shower-thermal" and "shower-thermal-thermal" quark recombination. These last three processes involving 
thermal quarks are illustrated in Fig. [5] 

The contributions from "shower-shower-shower" quark recombination, denoted as Dq SSS \ has the same expression 
as in vacuum [Ea. f5f)|> ] with the replacement of the triple quark distribution functions F qiq2q3 (zi, Z2, 23) in vacuum by 
its counterpart in medium F qiq2q3 (zi, Z2, Z3). Working along the same approach, one can also find out the contribution 
from "thcrmal-thermal-thcrmal" quark recombination, 



f\B(TTT) 1 n+) f d 2 qi±dz qi f d 2 q 2 ±dz q2 f 
[ B ' P } J 2(2^)3 J 2(2^)3 J 



ZB 



2(2tt) s 

ZB 



d 2 q 3 _ L dz q 
2(2tt) 3 



■fqi {qi±,z qiP + )f q2 (q 2 ± 

? Zq 2 P 



~qi z qi z q3 • 

<PB(qi±, — ;<Z2_L, — ;?3_L, 



ZB 



ZB 



-1 Z 

(2ir) 3 5( 2 \pB±-qi-L-q2±-q3±)S(zB-z qi -z q2 -z q3 ) , (104) 



where S q (zs) is the same factor associated with phase-space integration as in Eq. 197(1 . After extracting the spatial 
volume from the (5-function squared, one has 



V> 4 



d 2 p B \ 
(2tt) 3 



dz. 



2z B Jo 



2,1 dz, 



</2 



d 2 qi± f d 2 q 2 ± 



2z B J (2tt)3 J (2tt) 



fqi(qi±,z qi p + )f q2 (q2±,z q2 p + )f q3 (pB±-qij_-q2±,(zB-z qi -z q2 )p + ) 



ZB Z B 



), 



(105) 



which corresponds to an invariant inclusive baryon spectrum from thermal quark recombination, 



(2tt) ; 



dN B ^ TTT ) 
dptd 2 p B ± 



V / dx 



qx 



dx 



'12 



d 2 qu 



d 2 q 2 A 



2(2tt) 3 J 2(2tt) 3 



;fqAqi±,x qiP +)f q2 (q 2 ± 



x fq 3 (PBX-qi±-q2±, (l-x qi -x q2 )pg)R B {qi±, x qi ; q 2 ±, x q2 ) 



(106) 



The derivation of the contribution from the "shower-thermal-thermal" quark recombination to the baryon fragmen- 
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tation function is also straightforward. One has, 



D B ^(z B ,p + ) = J2 
93 gb 



d 2 qi± f d 2 q 2 ± 



(2tt) 3 
1 



(2tt) 3 



dz„ 



2z B 



Zqi dz 



2^/«l fe-L' 2 9lP + )/92(<?2_L, Z 9 2-P + ) 



(i - - z q2 /z B y 



z z — 

tRb(Qi±, —\<12±, —)Fa 3 (zB~z gi ~z q2 ), 



ZB 



ZB 



(107) 



where the summation is over different quark flavors of q% G B among the constituent quarks of the baryon B. 

In deriving the contribution from "shower-shower-thermal" quark recombination, one encounters again the problem 
of interference between amplitudes of different shower quark recombination. We have to make the same approximation 
that one can neglect the interference term and the final result is proportional to the diagonal contribution with a 
constant Cb2- Therefore, we have 



D^ sst Hzb,P + ) 



(9l92)SS 

1 



d 4 p 



8 z B - 



pi 



d 4 xe- ipx 
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(9l92)GB 
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(2tt)3 Ax kl J (2tt)3 2 



— f Q3 (q±,x q p'j 3 )\ip B (q±,x q ;k 1 ±,x kl - 1 k 2 ±,x k2 
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Tr 



77T 2^ k ( (^(0) | S; fci ± , a: fcl ; fc 2 _L , i£fc 2 ) (fc 2 _L , 

IV B S 



x k2 ;kij_,x kl ;S\ip(x)))k 



(108) 



where k 2 ± = pB±-qi±-ki±, k 2± = PB±-qi±-k' 1± , x k2 = l—x gi —Xk 1 , x' k2 = l—x gi —x' kl , Cb2 is a constant with the 
dimension of momentum and the summation is over all possible (qiq 2 ) £ B quark pairs among the three constituent 
quarks of the baryon B. Changing the variables x kl ,x k2 ,x qi to z kll z k2 , z qi and neglecting the intrinsic momentum 



in the baryon wavefunction, \lpb(q±, x q ] ki±, x kl ; fc 2 j_, x k2 



R B (qj_,x q ;0±,x kl ), we have 
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.fq 3 (<l^iZqP + ) 



R B (q x ,^;0 x ,—)F^(z,z B 



ZB 
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(109) 



The diquark distribution F^ l9 "(zi, z 2 ) is the same as in the "shower-shower" contribution to the meson fragmentation 
in Eq. Q29[). which is similarly defined as in the vacuum except that the thermal average in principle should include 
medium effects such as induced radiation and absorption. 

Finally, the effective baryon fragmentation function is the sum over the contributions resulting from "shower-shower- 
shower" , "shower-shower-thermal" , "shower-thermal-thermal" and "thermal-thermal-thermal" constituent quark re- 
combination, 



Df(z B ,P + ) = D^ SSS \z B ) +Df< SST >(z B ,P + ) +D^ STT \z B ,p + ) +D^ TTT \z B ,p + ) , 



(110) 



though the "thermal-thermal-thermal" contribution is not correlated with the initial parton jet and thus is only the 
background from the hadronization of the thermal medium. 



VI. CONCLUSION recombination during the hadronization of the parton jet 

together with a quark-gluon plasma. 

In this paper, we have studied medium modification to We started with a formulation of the vacuum par- 
the effective parton fragmentation functions due to quark 
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ton fragmentation functions in the parton operator form 
within a constituent quark model. In this model, con- 
stituent quarks are the effective degrees of freedom during 
the hadronization of both the parton jet and the quark- 
gluon plasma, which is essentially a thermalized gas of 
constituent quarks. Final hadrons are composed of con- 
stituent quarks with given wavefunctions. We showed 
that for sharply peaked hadron wavefunctions, one can 
neglect the interference between amplitudes of hadron 
formation from recombination of constituent quarks with 
different momenta. Consequently, we were able to cast 
the meson (baryon) fragmentation functions into a con- 
volution of the recombination probability and constituent 
diquark (triquark) distribution functions of the fragment- 
ing parton. The recombination probability is determined 
by the hadron's wavefunction in the constituent quark 
model. The diquark or triquark quark distribution func- 
tions of the fragmenting parton are defined as the over- 
lapping matrices between parton field operator and the fi- 
nal constituent quark states, similarly as the dihadron or 
trihadron fragmentation functions defined as the overlap- 
ping matrices between parton operator and final hadron 
states. We derived the QCD DGLAP evolution equa- 
tions for the diquark and triquark distributions functions 
which are a little different from that of dihadron and tri- 
hadron fragmentation functions because of the kinematic 
constraints imposed by the hadrons' wavefunctions in the 
constituent quark model. We further discussed possible 
connections between triquark, diquark and single quark 
distribution functions through sum rules. 

Working within the framework of field theory at finite 
temperature, we extended the formulation of parton frag- 
mentation functions in terms of quark recombination to 
the case of jet fragmentation in medium, assuming that 
the medium is hadronizing together with the parton jets. 
Replacing the vacuum expectation in the overlapping ma- 
trices by the thermal average, we were able to derive the 
medium modification to the fragmentation functions via 
quark recombination. The medium modification comes 
not only as the medium modification of the diquark or 
triquark distribution functions in the recombination of 
constituent quarks from the fragmenting parton, called 
shower quark recombination, but also as additional con- 



tributions from recombination between shower and ther- 
mal constituent quarks. We also obtained contributions 
from recombination of thermal quarks within the same 
formalism. However, such contributions of thermal quark 
recombination are disconnected with parton jet fragmen- 
tation and therefore are just a thermal background to jet 
fragmentation processes in the medium. 

The formulation and derivation of jet fragmentation 
functions in terms of quark recombination in the vac- 
uum do not simplify the theoretical description of the jet 
fragmentation processes. It is only a model of the non- 
perturbative process of hadronization. However, once ex- 
tended to jet fragmentation in medium, it provides a nat- 
ural framework for the description of jet hadronization 
in medium and the study of effective medium modifica- 
tion to the fragmentation functions. Within this formal- 
ism, one can not only include contributions to jet frag- 
mentation from the recombination of constituent quarks 
from the parton jet and the thermal medium but also the 
medium modification of the quark distributions of the jet 
parton due to induced gluon bremsstrahlung and absorp- 
tion. This framework will effectively unify parton energy 
loss via induced gluon radiation and quark recombination 
processes during hadronization. 
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